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The transformation of monolayer islands into bilayer islands as a first step of the overall two- 
dimensional to three-dimensional (2D-3D) transformation in the coherent Stranski-Krastanov mode 
of growth is studied for the cases of expanded and compressed overlayers. Compressed overlayers 
display a nucleation-like behavior: the energy accompanying the transformation process displays a 
maximum at some critical number of atoms, which is small for large enough values of the misfit, and 
then decreases gradually down to the completion of the transformation, non-monotonically due to 
the atomistics of the process. On the contrary, the energy change in expanded overlayers increases up 
to close to the completion of the transformation and then abruptly collapses with the disappearance 
of the monoatomic steps to produce low-energy facets. This kind of transformation takes place only 
in materials with strong interatomic bonding. Softer materials under tensile stress are expected 
to grow predominantly with a planar morphology until misfit dislocations are introduced, or to 
transform into 3D islands by a different mechanism. It is concluded that the coherent Stranski- 
Krastanov growth in expanded overlayers is much less probable than in compressed ones for kinetic 
reasons. 

PACS numbers: 68.35.Md, 68.35.Np, 68.65.Hb, 68.43.Hn 



I. INTRODUCTION 

Understanding the mechanism of transition from a pla- 
nar two-dimensional (2D) thin film to a three-dimensio- 
nal (3D) island morphology in the heteroepitaxy of highly 
mismatched materials is of crucial importance for the 
growth of self-assembled quantum dots in nanoscale tech- 
nology. In the coherent Stranski-Krastanov (SK) mode 
of growth, dislocation-free 3D islands develop on top of 
a 2D wetting layer in order to relieve the misfit strain at 
the expense of an increase in surface energyi This mech- 
anism of strain relaxation is established in a multitude of 
systems of technological importance for the manufactur- 
ing of optoelectronic devices.^ Despite the huge amount 
of studies devoted to the evolution of the cluster shape, 
many aspects, in particular the very beginning of the 2D- 
3D transition, still remain unclear. 

The first theoretical concept for the transition from 
a 2D layer to faceted 3D islands included a nucleation 
mechanism as a result of the interplay between the sur- 
face energy and the relaxation of the strain energy rela- 
tive to the values of the wetting layer. Irreversible 3D 
growth was predicted to begin above a critical volume, 
overcoming an energetic barrier whose height is inversely 
proportional to the forth power of the lattice misfit. 

Mo et al. observed Ge islands representing elongated 
pyramids ("huts") bounded by {105} facets inclined by 
11.3° to the substrate.'* These clusters were thought to be 
a step in the pathway to the formation of larger islands 
with steeper side walls ("domes" or "barns" )i^ Chen et 
al. studied the earliest stages of Ge islanding and found 



that Ge islands smaller than the hut clusters do not in- 
volve discrete {105} facets.^ This result was later con- 
firmed by Vailionis et al. who observed the formation of 
3-4 monolayers-high "prepyramids" with rounded bases 
which existed over a narrow range of Ge coverages in 
the beginning of the 2D-3D transformation>i Sutter and 
Lagalljf^ assumed that faceted, low-misfit alloy 3D is- 
lands can result from morphological instabilities (ripples) 
that are inherent to strained filmS)2ii2iii thus suggesting 
that 3D islands can be formed without the necessity to 
overcome a nucleation barrier. Similar views were si- 
multaneously expressed by Tromp et ali^ Tersoff et al. 
developed further this idea suggesting that the transi- 
tion from the initial smooth "prepyramids" to faceted 
pyramids can be explained by assuming that the polar 
diagram of SiGe alloy islands allows the existence of all 
orientations vicinal to (001) with the first facet being 
{105}.^ 

In order to explain the experimental observation that 
SiGe alloy films roughen only under compressive stresses 
larger than a critical value of 1.4 %, Xie et al. assumed 
that the smallest 3D islands have stepped rather than 
faceted side surfacesii^ They noted that the steps on the 
SiGe(OOl) vicinals are under tensile stress and their ener- 
gy of formation is lowered by the compressive misfit but 
increased by the tensile strain [see also the discussion in 
Refs. fill) and (16)]. As a result, the step formation and 
in turn the roughening are favored by the compressive 
misfit. It is worth noting that a barrierless evolution 
of stepped islands was predicted by Sutter and Lagally 
under the assumption that the slope of the side walls of 
the stepped islands increases continuously from zero to 
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11.3°.^ 

Priester and Lannoo suggested that 2D islands of 
monolayer height appear as precursors of the 3D is- 
lands. In addition, it was established that the mini- 
mum-energy pathway of the 2D-3D transition has to con- 
sist of a series of intermediate states with thicknesses 
increasing in monolayer steps and which are stable in 
separate intervals of volume. The first step in this trans- 
formation should be the rearrangement of monolayer into 
bilayer islandss^^ Khor and Das Sarma found by Monte 
Carlo simulations that during the rearrangement, the 
material for the bilayer island comes almost completely 
from the original monolayer island, the bulk of the ma- 
terial for the three-layer island comes from the original 
two-layer island, etci^ 

Moison et al. reported that the coverage suddenly de- 
creases from about 1.75 ML to 1.2 ML when 3D InAs 
islands begin to form on GaAs.^^ The same phenomenon 
has been observed by Shklyaev et al. in the case of 
Ge/Si(lll)f22i These observations suggest a process of re- 
arrangement as mentioned above. Voigtlander and Zin- 
ner noted that 3D Ge islands have been observed in 
the same locations on the Si(lll) surface where 2D is- 
lands locally exceeded the critical thickness of the wet- 
ting layer of two bilayers.^^ One-monolayer thick InAs 
islands were suggested to act as precursors for formation 
of thicker structures on GaAs.^"' The simultaneous pre- 
sence of stable one, two, three or four monolayers-thick 
islands has been observed in heteroepitaxy of InAs on 
InP and GaAsM*21 

In this paper we studied the earliest stages of growth 
of thin films in the coherent (dislocation-free) Stranski- 
Krastanov mode. We considered the instability of the 
planar growth against clustering by focussing on the con- 
servative (i.e., without considering further deposition) 
mono- to bilayer transformation as a first step of the over- 
all 2D-3D transition, or the beginning of the formation of 
the "prepyramids" mentioned above. We found that this 
transformation is a true nucleation process in compressed 
overlayers, in the sense that a critical nucleus of the sec- 
ond layer is initially formed and then grows further up to 
the complete mono-bilayer transformation. The energy 
associated with the transformation thus reaches a maxi- 
mum and then starts a decreasing trend. This is not the 
case in expanded overlayers, where the energy tends to 
increase up to very close to the completion of the trans- 
formation and then steeply decreases at the very end. 
The main result of this study is that coherent Stranski- 
Krastanov growth in expanded overlayers is much less 
probable than in compressed ones. 



II. MODEL 

We consider an atomistic model in 2 -I- 1 dimensions. 
The 3D crystallites have /cc structure and (100) surface 
orientation, thus possessing the shape of truncated squ- 
are pyramids. We found that monolayer- high elongated 



islands always have higher energy than square islands 
with the same number of atoms, as expected from the 
symmetry of the lattice and the isotropy of the interac- 
tion potentials (see below) . The lattice misfit is the same 
in both orthogonal directions. We consider interactions 
only in the first coordination sphere; inclusion of further 
coordination spheres does not alter qualitatively the re- 
sults. 

The choice of crystal lattice and interaction potential is 
more appropriate for the heteroepitaxy of (close-packed) 
metals on metals rather than for that of semiconductor 
materials. As a consequence, properties that depend cru- 
cially on the strong directional bonding characteristic of 
semiconductors cannot be addressed by our model. Some 
examples are: the dependence of the shape of GeSi/Si 
dots on volum o'^i^i^ as discussed above; the observation 
of lens-shaped22i^ and pyramidal^i dots and even of co- 
existence of both types^S. in InAs/GaAs, the other well- 
studied system (for a recent review see Ref. i33) or the 
cases where the accommodation of the lattice misfit of a 
given material on different crystallographic faces of the 
same substrate takes place by other mechanisms (also 
found for InAs/GaAs) ^24 where additional aspects as the 
presence of different surface reconstructions affect the 
thermodynamical balance of surface energies as well as 
the diffusion kinetics and, as a consequence, the nucle- 
ation behaviour and the growth mode. As our aim is to 
study the "reversible" minimum-energy pathway of the 
transition from metastable states to the ground state of 
a given system, the exact particularities of the model are 
not likely to play a crucial role and we expect the same 
qualitative behavior for any crystal lattice, crystal shape 
and interatomic potential. 

We have performed atomistic calculations making use 
of a simple minimization procedure. The atoms inter- 
act through a pair potential whose anharmonicity can 
be varied by adjusting two constants /i and v > v) 
that govern separately the repulsive and the attractive 
branches, respectively 



V{r) = Vo 



(1) 



where b is the equilibrium atom separation. For ^ — 2v 
the potential ^ turns into the familiar Morse form. 

A static relaxation of the system is performed by allow- 
ing each atom to displace in the direction of the force, i.e., 
the gradient of the energy with respect to the atomic co- 
ordinates, in an iterative procedure until all the forces fall 
below some negligible cutoff value. As we were only in- 
terested in the 2D-3D transformation of isolated islands, 
the calculations were performed under the assumption 
that the substrate (the wetting layer) is rigid; the atoms 
there are separated by a distance a. The lattice misfit is 
thus given by e = (6 — a)/a. 

For the study of the mechanism of mono-bilayer trans- 
formation, we assume the following imaginary model 
processs2J, atoms detach from the edges of monolayer is- 
lands, which are larger than the critical size for the mono- 
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FIG. 1: Schematic process for the evaluation of the activation 
energy of the monolayer-bilayer transformation. The initial 
state is a square monolayer island. The intermediate state is 
a monolayer island short of some number of atoms which are 
detached from the edges and placed in the second level. The 
final state is a truncated bilayer pyramid. 



to bilayer transformation and thus unstable against 
bilayer islands, diffuse on top of them, aggregate and give 
rise to second-layer nuclei. These grow at the expense of 
the atoms detached from the edges of the lower islands. 
The process continues up to the moment when the up- 
per island completely covers the lower-level island. To 
simulate this process, we assume an initial square mono- 
layer island, detach atoms one by one from its edges and 
locate them on top and at the center of the ML island, 
building there structures as compact as possible (Fig. ^1. 
The energy change associated with the process of trans- 
formation at a particular stage is given by the difference 
between the energy of the incomplete bilayer island and 
that of the initial monolayer island. This is in fact a 
conservative version of the mechanism observed by Khor 
and Das Sarma in dimensions^ 



III. RESULTS 



A. Stability of Monolayer Islands 
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FIG. 2: Total energy per atom of mono- and bilayer islands 
at (a) positive {+4.0 %) and (b) negative (-11.0 %) values of 
the misfit as a function of the total number of atoms. The 
potential is of the form given by eq. with fi — 16 and 
u = 14. 



In our previous work in (1+1)D models, it was estab- 
lished that monolayer-high islands are stable against bi- 
layer islands up to some critical volume or number of 
atoms N12; in turn, bilayer islands are stable against tri- 
laycr islands up to another critical number > .^12; 
g-j-^, 18,19 rpj^g mono-bilayer transformation was considered 
as the first step of the overall 2D-3D transformation and a 
critical misfit was determined from the misfit dependence 
of Ni2- The latter was found to increase with decreas- 
ing misfit diverging at a critical value £12. Above the 
critical misfit, the coherent Stranski-Krastanov mode is 
favored against the layer-by-layer growth followed by in- 
troduction of misfit dislocations. The opposite is true 
below the critical value. Whereas this critical behavior 
is clearly pronounced with compressive strain, it is much 
smoother in expanded overlayers. It is worth noting that 
the existence of critical misfit was observed in a series of 
heteroepitaxial systems ii^i2S*^2i2S 

In the present work, using more realistic (2-t-l)D mod- 
els, we found a larger difference in the behavior of ex- 
panded and compressed overlayers. Figure |21 shows that 
the total energies of mono- and bilayer islands under ten- 



sile stress containing the same total number of atoms are 
very close to each other compared with the corresponding 
behavior in compressed overlayers. As will be shown be- 
low this leads to the conclusion that even for TV ^ N12, 
determined as the crossing point of the curves corre- 
sponding to monolayer and bilayer islands in Fig. |21 the 
probability of the 2D-3D transformation remains nearly 
equal to the probability of the reverse 3D-2D transfor- 
mation. 

It turns out that the misfit dependence of is very 
sensitive to the value of the force constant 7 = /iz^Vb of 
the interatomic bonds, particularly in expanded overlay- 
ers (Fig. Decreasing fi and v {Vq is assumed equal 
to unity) in such a way that the ratio ^/v is kept con- 
stant (in this case equal to 8/7), shifts the intersection 
points N12 to larger absolute values of the misfit. As a 
result, a critical size 7Vi2 in compressed overlayers exists 
practically for all values of 7 whereas in overlayers under 
tensile stress, A^i2 shifts to so large values of the misfit 
that they effectively disappear below some critical value 
of 7. 
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FIG. 3: Critical island size N12 (number of atoms) as a func- 
tion of the lattice mismatch at different values of the force 
constant 7 — fifVo- Potentials of the form given by eq. Q 
were used, with fi/i/ = 8/7 (Vb is taken equal to unity). As 
seen, coherent 3D islanding is favored in expanded overlayers 
only in "stiff er" materials. 



It was established that in the case of a force constant 
of an intermediate value (fj, = 2u = 12), N12 disappears 
(the monolayer islands are always stable against bilayer 
islands), but N13, N14, N23--- still exist. This points to a 
novel mechanism of 2D-3D transformation which differs 
from the consecutive formation of bilayer, trilayer, etc. 
islands, each from the previous one. The new mecha- 
nism obviously consists of the formation and 2D growth 
of bilayer islands on top of the monolayer island, thus 
transforming the initial monolayer island directly into a 
trilayer island. At even smaller values of 7, the critical 
values 7Vi3,7Vi4 etc. consecutively disappear, suggest- 
ing a generalized mechanism of 2D-3D transformation in 
which the monolayer islands transform into thicker is- 
lands. This multilayer 2D mechanism will be a subject 
of a separate study. In this paper we focus on the layer- 
by-layer 2D-3D transformation. 

We considered also the stability of monolayer islands 
against bilayer islands with a different slope of the side 
walls. It was found that N12 is smaller if the slope of 
the side walls is the steepest one (60°) for this lattice in 
comparison with flatter islands (Fig. 01 . This is due to 
the fact that in crystals with steeper side walls, the strain 
relaxation is more efficient than in flatter islands. This is 
in contradiction with the experiments in semiconductor 
growth in which islands with side walls of a smaller slope 
than that of the first facet are initially observedALi^ 
In any case this means that we can exclude the flatter 
islands from our consideration. 

We conclude that for some reasonable degree of an- 
harmonicity (e.g. /i = 2i/ in our model), monolayer is- 
lands become unstable against bilayer islands thus mak- 
ing possible the 2D-3D transformation by the layer-by- 
layer mechanism only at strong enough interatomic bond- 
ing. Soft materials are expected to grow either with a pla- 



FIG. 4: Misfit dependence of the critical sizes of mono- 
and bilayer islands with different shapes, given by different 
angles of the side walls: 60° and 30° {fi — 21/ = 12). The 
critical misfit £12 is shown by the vertical dashed line. 



nar morphology until misfit dislocations are introduced 
or to transform into 3D islands by a different, multilayer 
2D mechanism. 



B. Mechanism of 2D-3D transformation 

FigureElshows typical transformation curves of the en- 
ergy change as a function of the number of atoms in the 
upper island for (a) positive and (b) negative misfits. It 
is immediately seen that in a compressed over layer, the 
transformation curve for AG has the typical shape of a 
nucleus formation curve: it displays a maximum AG^ax 
for a cluster consisting of a small number of atoms rimax 
and then decreases beyond this size up to the completion 
of the transformation. The atomistics of the transfer pro- 
cess (i.e. the completion of rows in the upper level and 
their depletion in the lower one) is responsible for the 
jumps (the non-monotonic behaviour) in the curve. 

In the case of expanded overlayers, the energy change 
increases up to a large number of atoms (again non- 
monotonically due to the atomistics) and then abruptly 
decreases at the end of the transformation. No true max- 
imum is displayed. The energy change becomes negative 
only after the transfer of the last several atoms. Com- 
paring the largest value of the energy with the energy 
at the transformation completion leads to the conclusion 
that the probabilities of the direct and reverse transfor- 
mations are nearly equal. 

Figure depicts the evolution of the height of the bar- 
rier AGmax as a function of misfit (in expanded over- 
layers, this is the highest value reached before the col- 
lapse of the energy). The figures at each point show the 
number of atoms in the cluster at the maximum of the 
transformation curve. As seen, in the case of compressed 
overlayers, AGmax decreases steeply with increasing mis- 
fit in a way similar to the decrease of the work required 
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FIG. 5: Transformation curves representing the energy 
change in units of the bond energy Vb as a function of the 
number of atoms in the upper level for (a) positive (+2.5 %) 
and (b) negative (-7.0 %) values of the misfit. The number of 
atoms in the initial monolayer island (A'^o = 841 = 29 x 29) 
is chosen in such a way that the resulting truncated bilayer 
pyramid is complete (21 x 21 = 441 atoms in the lower and 
20 X 20 = 400 atoms in the upper level); ^i^2u = 36. 



for nucleus formation with increasing supersaturation in 
the classical theory of nucleation.^* Assuming a depen- 
dence of the kind AG = Xe~" where K \s a. constant 
proportional to the Young modulus (or the force con- 
stant 7) and £ is the lattice misfit, we found n — 4.29 
for 11 = 12,^ = 6, and n = 4.75 for fi — 36,1^ — 18. 
It is worth noting that assuming 3D nucleation on top 
of the wetting layer, Grabow and Gilmer predicted a 
value n = 4 for small misfits (large nuclei) assuming 
that AGraax is inversely proportional to the square of 
the supersaturation, which in turn is proportional to the 
square of the lattice misfit 4^ Note that the same expo- 
nent of four was obtained also by Tersoff and LeGouesi^ 
Obviously, in our case the exponent rt is a complicated 
function of the force constant in the interatomic bonds 
but the value of the exponent is of the same order. 

The misfit behavior of the critical nucleus size Umax 
is also similar to that found in the theory of nucleation. 
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FIG. 6: Height of the energetic barriers in units of Vb as a 
function of the absolute value of the lattice misfit. The figures 
at each point show the number of atoms in the critical nucleus 
n,nax- The initial island size was 29x29 = 841 atoms. The 
round symbols were calculated for jj. = 2v = 36, the squares 
for fi^2u = 12. 
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FIG. 7: Misfit dependence of the critical island size A'^12, 
the critical nucleus size Umax, (both expressed in number of 
atoms), and the nucleation barrier height AGmax (in units 
of Vb) for compressed overlayers and ^ — 2v = 12. The last 
two magnitudes were computed for islands of an initial size of 
20x20 = 400 atoms. 



The nucleus size decreases with increasing misfit, reach- 
ing eventually only one atom just as the nucleus size on 
the supersaturation in nucleation theoryjSS, It is inter- 
esting to note that the critical nuclei contain always a 
number of atoms which exceeds by one atom the size 
of a compact cluster, and is not one atom short of be- 
ing a compact cluster, as might be expected from too- 
simplistic energetic considerations on the basis of bond- 
counting arguments. The reason is that the highest de- 
parture from compactness giving the maximum of the 
AG curve is achieved when the additional atom creates 
two kink positions for the attachment of the next row of 
atoms. 
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FIG. 8: Variation of the in-plane strain energy and out- 
-of-plane interaction energies during the mono-bilayer trans- 
formation process in (a) compressed (e = -1-2.5 %) and (b) 
expanded (e = —7.0 %) overlayers with n = 2i> = 36. The 
initial sizes of the islands were 29 x 29=841 atoms. The inserts 
show the curves at the beginning of the transformation with 
enlarged scales. 



On the contrary, in expanded overlayers, the number 
of atoms at which the transformation curve reaches its 
highest value does not depend on the misfit (Fig.El), de- 
monstrating the non-nucleation behavior of the process. 
This number is roughly equal to the number of atoms 
which completes the upper level minus the number of 
atoms required to build the last four edge rows of the 
upper level in order to produce four facets. A special 
feature of the above results is that the barriers for ex- 
panded overlayers are in general larger than those for 
compressed overlayers. Having in mind that the typi- 
cal time needed for the transformation to occur is in- 
versely proportional to exp{—AGmax/kT), we have to 
expect much longer times for the occurrence of the 2D- 
3D transformation in expanded overlayers as compared 
with compressed ones particularly at larger values of the 
misfit. Thus limitations of kinetic origin (astronomically 
long times for second layer nucleation) are expected in 



expanded overlayers and at small misfits in compressed 
overlayers. 

Figure [T] compares the misfit dependences of the crit- 
ical island size N12, the critical 2nd-layer nucleus size 
TT-max and the height of the nucleation barrier AGmax for 
compressed overlayers. The tree curves display a similar 
behavior increasing steeply with decreasing misfit, with 
N12 showing a critical behavior at the critical misfit £12- 
As will be discussed below, the 2D-3D transformation 
will be inhibited for kinetic reasons at values of the mis- 
fit not sufficiently larger than £12. 



IV. DISCUSSION 

In compressed overlayers the atoms interact with their 
in-plane neighbours through the steeper repulsive branch 
of the interatomic potential. This means that compressed 
overlayers are effectively "stiffer" than expanded over- 
layers. Then compressed pseudomorphic overlayers will 
contain a larger amount of elastic energy than expanded 
pseudomorphic ones of the same thickness and bonding 
strength. Another consequence is that the accumulation 
of strain energy with thickness in compressed overlay- 
ers will be steeper than in expanded overlayers. Or, the 
strain energy of a bilayer will differ considerably from 
that of a monolayer from both sides of the critical size iVi2 
compared with expanded islands as seen in Fig. [3 This is 
the reason why expanded overlayers require greater force 
constants in order for the monolayer islands to become 
unstable against bilayer islands. This is also the reason 
why the critical sizes A'^12 are larger in expanded than in 
compressed overlayers of the same force constant as seen 
in Fig.13 

Figure |31 illustrates the central result of our study. It 
shows the energy of transformation of compressed and 
expanded monolayer islands. In compressed overlayers 
the overall fall of the energy begins when a small cluster 
(three atoms in this particular case) is formed in the sec- 
ond level. This is a typical phase transition of first order - 
nucleation followed by irreversible growth. In expanded 
overlayers the energy increases up to nearly the end of 
the transformation and then abruptly collapses. This 
collapse of the energy is connected with the transfer of 
the remaining last edge rows of atoms which leads to the 
coalescence of the lower and upper steps to produce four 
side facets. In the particular transfer process considered 
in the calculation, there are two regions, one before the 
beginning of the collapse, the second between the two 
sub-collapses, where the total energy rises slightly; this 
is due to the energetic cost of repulsion between steps 
that arc very close, separated only by a single atomic 
row. The reason for this "non-nucleation" behavior is 
the effectively weaker expanded bonds which results in 
relatively close energies of the monolayer and bilayer is- 
lands. With increasing size of the second level cluster, 
the misfit strain is not as effectively relaxed as in the 
case of compressed islands and the collapse of the energy 
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is due to the replacements of isolated repulsing steps by 
low-energy facets. 

The different transformation behavior can be under- 
stood on the basis of the above considerations account- 
ing in addition for the finite size of the islands. During 
the transformation of the monolayer islands we should 
expect a relaxation of the in-plane strain, which is the 
driving force for the transformation and an increase of 
the total energy of interaction between the lower level of 
the island and the substrate and between both levels of 
the islandiSiia 

The physics behind the process of mono-bilayer trans- 
formation is the same as behind the formation of 3D 
islands on the wetting layerii^ Formation and growth 
of new steps and relaxation of the strain stored in the 
monolayer island compete in the process. In addition, 
the steps in the first and second levels repel each other 
and the associated repulsion energy increases towards the 
completion of the transformation to disappear when the 
steps coalesce to give rise to microfacets. 

Epitaxial strain is relaxed at the islands edgesn^ Edge 
atoms are displaced from the bottoms of their respective 
potential troughs giving rise to relaxation of the bonds 
parallel to the surface (in-plane stress relaxation). The 
displaced atoms loose contact with the substrate atoms, 
which leads to an increase of the out-of-plane energy 
of interaction with the underlying substrate. ^"^-^^ Dur- 
ing the transformation, new edges on top of the initial 
monolayer island are formed and the total length of the 
edges increases as AL{N2) = 4(VA^o - N2 + \^-VN^) 
where A^o and N2 are the number of atoms in the initial 
monolayer island and the current number of transferred 
atoms in the second level, respectively. Note that the en- 
ergy of repulsion between the edges bounding the lower 
and upper islands is implicitly accounted for by both the 
in-plane and out-of-plane energies. Then a larger number 
of atoms are displaced from the bottom of their respec- 
tive potential troughs during the transformation process 
and the total in-plane strain-relaxation energy decreases. 
Simultaneously, the out-of-plane interaction energy in- 
creases. Owing to the weaker attractive forces in ex- 
panded overlayers, only a small number of bonds close to 
the edges are relaxediiS Most of the bonds at the center 
of the islands are strained to fit the underlying wetting 
layer. As a result, the average relaxation in expanded 
islands is smaller than in compressed islands, where even 
bonds at the center of medium-sized islands are partly 
relaxed. In compressed overlayers, the decrease of the 
in-plane strain energy rapidly overcompensates the in- 
crease of the out-of-plane interaction energy which re- 
sults in a nucleation-like transformation curve (Fig. (SJ. 
In expanded overlayers, the absolute value of the total in- 
plane strain energy is smaller than the out-of-plane inter- 
action energy with the exception of the final stage when 
the monolayer-high steps disappear to produce facets of 
small surface energy. 

The typical time required for the appearance of a 



second-layer nucleus is inversely proportional to the nu- 
cleation frequency lu — S12K exp[—AGmax/kT), where 
S12 = a?'Ni2 is the area of the critical monolayer island 
and K is the pre-exponential of the nucleation rate. As 
seen in Fig. |H| in the case of fi — 2v — \2, the barrier 
height increases approximately 5 times in an interval of e 
of 2.5 % whereas the number of atoms in the critical nu- 
cleus rimax increases nearly 70 times. For a greater force 
constant — 2v — 36), the increase of both AGmax and 
^max is even larger: 20 and 110 times, respectively, in 
a smaller misfit interval of about 1.5 %. The energy to 
break a first-neighbor bond, Vq, for most semiconductor 
materials is of the order of 2 to 2.5 eV (the enthalpy of 
evaporation is of the order of 4 to 5 eV) . Assuming N12 is 
of the order of 100 - 120 atoms we could expect a mono- 
bilayer transformation to take place at misfits for which 
AGmax/kT < 15 — 20 {n„iax < 3). The reason is that 
the pre-exponential K in 2D nucleation rate from va- 
por is usually of the order of lO^*' cm^^s~^i^ Otherwise, 
due to the exponential dependence, times of the order of 
centuries would be required for second-layer nucleation^^ 
Thus, although in compressed overlayers second-layer nu- 
cleation can be expected for thermodynamic reasons at 
misfits above £12, a real 2D-3D transition can only take 
place at even larger misfits or higher temperatures for 
kinetic reasons. As the height of the transformation bar- 
riers in expanded overlayers is always greater than sev- 
eral times Vb, the mono-bilayer transformation should be 
strongly inhibited for kinetic reasons. 

We conclude that the case of a layer-by-layer mecha- 
nism for the 2D-3D transformation is expected only in 
compressed overlayers at misfits sufficiently larger than 
£12. The reason is that the mechanism of the mono- 
bilayer transformation is nucleation-like due to the inter- 
play of relaxation of the in-plane strain, which is pro- 
portional to the total edge length and the increase of 
the total edge energy and repulsion between the edges. 
The transformation curve in expanded overlayers shows 
a "non-nucleation" behavior characterized by an overall 
increase of the energy up to the stage when the single 
steps coalesce to produce low-energy facets. The latter is 
accompanied by a collapse of the energy. The maximum 
energy is large and 2D-3D transformation of expanded 
overlayers is not expected for kinetic reasons even for 
materials with strong interatomic bonds. Softer materi- 
als are expected to grow with a planar morphology until 
misfit dislocations are introduced, or to transform into 
3D islands by a different mechanism. 
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